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1.  Introduction 


Sensor  systems  operating  at  about  95  GHz  are  being  investigated  for  use  in 
smart  mimitions  designed  to  defeat  armored  vehicles  [1].  These  systems 
require  specific  information  on  the  millimeter  wave  (MMW)*  background 
environments  that  can  be  encotmtered,  in  order  to  predict  their  target- 
^  detection  performance  in  a  real  environment.  This  information  can  be  ob¬ 

tained  from  actual  MMW  backgroimd  measurements,  or  from  a  validated 
MMW  clutter  model.  The  data  presented  in  this  report  were  taken  to  assist 
in  the  validation  of  a  millimeter  wave  background  model  for  the  Smart 
*  Weapons  Operability  Enhancement  (SWOE)  program  managed  by  the  U.S. 

Army's  Cold  Regions  Research  and  Engineering  Laboratory.  The  aim  of 
the  SWOE  program  is  to  produce  validated  models  for  simulating  real- 
world  environments  for  both  the  infrared  (IR)  and  MMW  regions  of  the 
spectrum.  The  Army  Research  Laboratory  (ARL)  was  not  directly  involved 
in  the  MMW  modeling  effort,  but  to  support  the  modeling  work,  we  devel¬ 
oped  empirical  models  for  the  95-GHz  circular  polarization  radar  clutter 
data  obtained  over  a  6-  to  7-week  time  period  at  Yuma,  AZ,  and  Grayling, 
MI.  We  discuss  the  radar,  calibration,  and  measurement  procedures  in  sec¬ 
tion  2  of  this  report,  and  the  measured  data  in  section  3.  In  section  4,  we 
describe  the  empirical  modeling  procedure  and  results,  and  we  present 
our  conclusions  in  section  5. 


% 


*An  acronym  list  is  provided  at  the  end  of  this  document  to  assist  the  reader. 
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Ground-Based  Radar  Description  and  Measurements 


2.1  Description  of  the  Radar  System 

The  radar  used  for  these  measurements  has  been  discussed  in  detail  previ¬ 
ously  [2,3].  Only  a  brief  description  of  the  radar  as  configured  for  this  test 
will  be  presented.  The  radar  was  designed  for  obtaining  95-GHz  target- 
and  background-signature  data,  with  the  radar  located  on  a  tower  and 
computer-controlled  from  a  remote  ground  station.  Figure  1  shows  the  ra¬ 
dar  and  the  computer-controlled  pedestal  set  up  for  the  field  test  at  Gray¬ 
ling.  A  video  camera  that  is  boresighted  with  the  radar  records  the  radar 
scene  through  a  glass  window  at  the  upper-right-hand  comer  of  the  metal¬ 
lic  cover.  The  radar  power  supply  and  temperature-control  air  conditioner 
are  also  shown  in  the  photograph.  The  radar  is  protected  from  rain  or 
snow  by  the  tent  during  operation  and  storage.  Data  recording  and  control 
of  the  radar  are  performed  at  a  ground  station  located  in  an  11-m-long 
semitrailer.  Data  from  the  target  platform  and  a  meteorological  station 
may  also  be  recorded  at  that  location. 

The  basic  characteristics  of  the  ground-based  radar  are  summarized  below: 

•  Frequency — ^95.0  GHz, 

•  Two-coordinate  amplitude-comparison  monopulse, 

•  Polarization  agile — RHCP/LHCP  or  V/H  (selectable), 

•  Coherent  and  wide  bandwidth, 

•  Frequency  steppable — 1  to  512  steps,  1-,  2-,  5-,  10-,  20-MHz  stepsize, 

•  Peak  power— 45  W, 

•  Pulse  width — 100  ns, 

•  Pulse  repetition  frequency — 10  kHz, 

•  Antenna  beamwidth  (one  way) — 0.3  to  5°  (1.5°  for  these  tests), 

•  Receiver  IF  bandwidth — ^30  MHz, 

•  Receiver  noise  figure — 10  dB, 

•  Dynamic  range  (amplitude) — 66  dB,  and 

•  Data  recording — 16  channels  for  each  pulse. 

The  radar  is  a  completely  polarimetric,  two-coordinate,  amplitude- 
comparison  monopulse  system.  It  is  a  fully  coherent,  dual  polarization, 
pulsed  radar  that  operates  in  a  frequency-stepped  mode,  between  95.00 
and  95.64  GHz.  For  the  first  two  tests,  we  operated  the  system  with  64 
10-MHz  steps  and  a  transmit  PRF  of  10  kHz;  for  the  second  Grayling  test, 
we  used  128  5-MHz  steps.  The  peak  power  of  the  all-solid-state  transmitter 
is  45  W;  the  pulse  width  is  100  ns;  and  we  swept  the  location  of  the  end  of 


Figure  1.  Ground- 
based  95-GHz 
monopulse  radar  at 
Camp  Grayling,  MI. 


the  15-m-long  range  gate  sequentially  from  105  to  255  m.  We  switched  the 
polarization  from  pulse-to-pulse,  and  obtained  both  linear  and  circular 
data.  We  placed  primary  emphasis  on  the  collection  of  circular  polarization 
data,  so  there  is  approximately  three  times  more  data  for  circular  polariza¬ 
tion  than  for  linear  polarization.  The  radar  used  a  15-cm-diameter  lens  an¬ 
tenna,  with  a  one-way  beamwidth  of  1.5°. 

The  radar  DAS  has  six  I/Q  detectors  for  the  monopulse  sum  and  differ¬ 
ence  signals.  The  relative  phase  of  the  two-sum  charmels  and  their  log  am¬ 
plitudes  are  also  processed  in  the  DAS.  For  each  transmitted  pulse,  the  sig¬ 
nals  received  in  each  of  the  16  channels  are  sampled,  multiplexed,  and 
digitized  by  two  12-bit  A/D  converters.  The  multiplexed  data  are  trans¬ 
mitted  from  the  radar  via  a  fiber-optic  data  link  to  the  remote  computer- 
control /recording  system  at  a  rate  of  240  KB/s.  (The  latter  is  a  PC  that  in¬ 
cludes  an  80486  33-MHz  processor,  8  MB  of  memory,  and  a  1-GB  disc 
drive.)  Data  are  magnetically  recorded  on  2-GB  8  mm  cassette  tapes. 

The  radar  is  mounted  on  a  pedestal  that  can  be  computer-controlled  lo- 
^  cally  or  from  a  remote  station  to  scan  the  radar,  or  to  orient  it  to  a  fixed  po¬ 

sition.  The  pedestal  is  an  elevation-over-azimuth  type,  with  0.01°  pointing 
accuracy  and  scanning  rates  between  0.1  and  10°/s.  In  addition,  a  co- 
V  boresighted  video  camera  mounted  on  the  radar  head  provides  a  video 

record  of  the  scanned  area. 

2.2  Radar  Calibration  Procedure 

The  radar  calibration  procedure  is  complicated  due  to  the  number  and 
types  of  data  channels.  The  entire  procedure  has  been  reported  in  detail 
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previously  (ref.  2  and  3);  therefore,  we  will  discuss  only  a  brief  outline  of 
the  calibration  procedure  specific  to  these  measurements.  We  performed 
preliminary  calibrations  in  the  Laboratory  before  taking  the  radar  to  the 
field-measurement  site  at  Grayling,  and  we  did  a  series  of  premeas¬ 
urement  checks  and  calibrations  at  the  field  test  site.  The  basic  features  of 
these  checks  and  calibrations  [4]  included  taking  proper  account  of  gain 
imbalances  in  I/Q  detectors,  and  nonlinear  responses  in  detectors  and 
amplifiers. 

We  performed  the  radar  calibrations  at  the  field-measurement  site  using 
the  reflectors  listed  in  the  calibration  target  array  of  table  1.  We  chose  the 
range  and  the  RCS  of  reflectors  3, 4,  and  5  so  that  the  received-signal  levels 
from  the  reflectors  were  approximately  equal  and  in  the  linear  range  of  the 
receivers  (about  6  dB  below  A/D  saturation).  We  used  reflectors  1  and  2  for 
determining  a  sign  in  the  PDM  calibration  described  below,  for  circular 
and  linear  polarization,  respectively. 

We  used  the  PDM  calibration  technique  [5]  to  compensate  for  the  effects  of 
imperfect  cross-polarization  isolation  in  the  radar.  This  calibration  method 
is  intended  to  remove  the  effects  of  the  radar  system  parameters  from  the 
data  so  that  the  target  polarization  scattering  matrix  can  be  determined 
without  degradation  due  to  the  system.  It  can  be  used  for  both  linear  and 
circular  polarization  by  measuring  the  five  reflectors  listed  in  table  1.  For 
this  report,  we  carried  out  the  PDM  calibration  technique  separately  for 
both  circular  and  linear  polarization. 

For  PDM  calibration  purposes,  depending  on  the  polarization,  we  selected 
four  of  the  five  reflectors  and  took  200  ramps  of  data.  (A  ramp  of  data  con¬ 
sists  of  64  pulses  separated  in  frequency  by  10  MHz  that  are  transmitted 
with  1  polarization,  and  64  pulses  transmitted  with  the  orthogonal  polar¬ 
ization  for  the  10-MHz  step-size  data,  or  128  pulses  transmitted  at  each  po¬ 
larization  for  the  5-MHz  step-size  data.)  The  first  step  in  processing  the 
PDM  data  is  to  perform  the  1/Q  corrections.  We  generated  the  I/Q  correc¬ 
tions  by  characterizing  each  I/Q  detector  at  four  different  amplitude  levels 
at  22.5°  phase  intervals,  from  0°  to  360°.  This  I/Q  correction  reduces  the  re¬ 
sidual  amplitude  error  due  to  I/Q  imbalance  to  40  dB  below  the  signal 
level;  the  maximum  required  phase  correction  was  3°.  We  averaged  the  200 
ramps  of  data  before  calculating  the  PDM  coefficients.  To  properly  average 
the  coherent  channels  in  the  PDM  data,  a  constant  phase  for  each  fre¬ 
quency  is  required  to  calculate  the  correct  average  amplitude  and  phase. 
We  applied  a  motion-compensating  algorithm  to  the  data  before  averaging 
to  minimize  any  phase  variations  due  to  motion  of  the  calibration 


Table  1.  Calibration 
reflector  array. 


Reflector 

no. 

Reflector 

type 

RCS 

(dBsm) 

Range 

(m) 

1 

Griddecd  trihedral 

14 

160 

2 

Dihedral  22.5° 

18 

190 

3 

Trihedral 

23 

218 

4 

Dihedral  0° 

26 

248 

5 

Dihedral  45° 

26 

280 
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reflectors  during  the  acquisition  of  200  ramps  of  data.  We  then  calculated 
the  PDM  correction  parameters  for  each  frequency  and  for  each  polariza¬ 
tion  from  the  averaged  data.  We  carried  out  this  calibration  procedure  be¬ 
fore  and  after  all  data  runs,  and  used  the  average  of  these  calibrations  for 
calibrating  the  data  for  the  first  Grayling  test.  The  calibration  data  were  ac¬ 
quired  as  part  of  the  scan  data  for  Yuma  and  the  second  Grayling  test. 

The  cross-polarization  isolation  of  the  radar  was  between  25  and  30  dB 
across  the  frequency  band,  and  the  improvement  in  the  cross-polarization 
isolation  using  the  PDM  calibration  technique  was  between  5  and  10  dB, 
depending  on  the  frequency.  The  improvement  due  to  the  PDM  calibration 
technique  is  limited  by  the  clutter  return  from  the  antenna  sidelobes,  as 
well  as  the  reflected  signal  from  the  comer  reflector  mounting  poles.  The 
resulting  cross-polarization  isolation  of  the  system  was  improved  to  about 
35  dB,  which  is  more  than  adequate  for  measuring  MMW  backgrounds, 
since  the  cross-polarization  ratio  for  natural  backgrounds  was  found  to  be 
less  than  10  dB. 

We  obtained  antenna  radiation  patterns  at  the  test  site  using  the  various 
reflectors.  The  approximate  two-way  sum  channel  beamwidth  for  the 
15-cm  lens  used  for  the  measurements  reported  here  was  1.0°,  and  the  two- 
way  sidelobe  levels  were  less  than  40  dB.  We  determined  the  radar 
monopulse  boresight  angle  for  each  polarization  channel,  using  the  appro¬ 
priate  reflector.  We  foimd  that  the  boresights  for  the  different  polarization 
combinations  varied  slightly,  with  most  of  the  variations  in  the  azimuth. 
The  overall  boresight  angle  difference  for  the  RR  and  RL  channels,  for  ex¬ 
ample,  was  about  0.1°,  with  most  of  the  difference  in  the  azimuth.  We  set 
the  optical  boresight  on  the  video  camera  to  coincide  with  the  average  of 
the  individual  boresights. 

Clutter  Background  Measurement  Procedure 

The  high-intensity  measurement  areas  at  Grayling  and  Yuma  were  outfit¬ 
ted  with  a  full  array  of  meteorological  equipment  for  the  SWOE  program 
by  other  laboratories  in  order  to  obtain  standard  meteorological  data,  as 
well  as  measurements  of  solar  radiance,  soil  moisture,  soil  temperature, 
and  cloud  conditions.  In  addition,  snow  depth  and  snow  density  were 
measured  in  Grayling,  when  appropriate.  The  measurement  area  in  Gray¬ 
ling  was  100  by  100  m,  and  was  located  at  the  bottom  of  a  small  valley.  The 
radar  was  positioned  on  top  of  a  hill  overlooking  the  valley  so  that  the 
measurement  area  could  be  scanned  effectively.  For  Yuma,  the  radar  was 
mounted  on  top  of  a  12-ft-high  platform  overlooking  the  two  measure¬ 
ment  areas,  defined  as  the  Thermal  Imaging  Processing  System  (TIPS)  and 
Atmospheric  Sciences  Laboratory  (ASL)  fields  of  view.  These  measure¬ 
ment  areas  were  located  on  either  side  of  the  platform;  the  calibration  re¬ 
flectors  were  located  midway  between  the  two  measurement  areas. 

At  Grayling  1  and  Yuma,  we  scanned  the  measurement  areas  using  a  raster 
pattern  covering  an  azimuth  increment  of  1.2°,  an  elevation  increment  of 
1.0°,  and  a  range  interval  of  150  m  in  15  m  increments.  The  raster  pattern 
slightly  overscanned  the  areas  in  azimuth.  For  the  Grayling  II  test,  we  in¬ 
creased  the  azimuth  overscan  by  using  1.0°  azimuth  steps;  we  also  reduced 


the  elevation  steps  slightly.  The  total  azimuth  scan  varied  from  30°  to  50°, 
depending  on  the  test  site.  We  used  range  overlaps  to  keep  reductions  in 
the  measured  RCS  caused  by  the  sum  channel  beam  shape  function  to  a 
minimum  over  the  entire  area.  During  the  first  two  tests,  we  recorded  32 
ramps  of  data  for  each  azimuth,  elevation,  and  range  gate.  (As  discussed 
previously,  each  ramp  consists  of  128  transmitted  pulses.)  During  the 
second  Grayling  test,  we  reduced  the  frequency  step  size  to  5  MHz,  and 
transmitted  256  pulses  during  each  ramp.  This  step  allowed  us  to  use  anti¬ 
aliasing  techniques  to  generate  high-resolution  RCS  maps.  In  addition,  we 
recorded  only  8  ramps  of  data  for  each  cell  location.  At  each  angular  posi¬ 
tion,  we  swept  the  radar  range  gate  through  the  range  interval,  and  then 
incremented  the  azimuth  angle.  When  all  azimuth  angles  had  been  meas¬ 
ured  at  one  elevation  angle,  we  incremented  the  elevation  position.  The 
depression  angles  varied  from  5°  to  1°  for  both  tests.  We  continued  this 
process  until  we  had  measured  the  entire  scan  pattern. 

We  took  the  MMW  measurements  at  Grayling  I  over  25  2-hour  time  peri¬ 
ods,  randomly  selected  from  19  to  23  September  and  17  to  23  October.  Dur¬ 
ing  the  two  hours,  we  made  three  complete  scans  (or  data  runs)  with  circu¬ 
lar  polarization,  and  one  data  run  with  linear  polarization.  The  raw  data 
file  included  the  six  coherent  1/ Q  monopulse  data  channels,  the  two  loga¬ 
rithmic  sum  channels,  and  the  relative  polarimetric  phase  I/Q  channel. 
The  final  calibrated  data  set  includes  72  circular  and  22  linear  polarization 
runs.  The  Yuma  data  set  includes  201  circular  and  27  linear  runs  taken  over 
a  47-day  period,  while  the  Grayling  II  data  set  includes  268  circular  and 
140  linear  runs  taken  over  a  43-day  period.  Since  the  same  cells  were  meas¬ 
ured  again  and  again  for  different  environmental  conditions  and  times,  the 
data  provide  a  good  basis  for  validation  of  a  SWOE  MMW  background 
model. 


3.  Data  Discussion 


3.1  RCS  Variations 

In  order  to  evaluate  the  effects  of  various  environmental  conditions  on  the 
measured  RCS  for  the  three  different  field  exercises,  we  selected  a  small 
subset  of  the  measured  cells  from  the  radar  scan  pattern  from  each  test  for 
further  analysis,  and  identified  them  so  that  each  background  type  was 
represented.  The  background  types  in  the  selected  cells  were  grass  and 
bushes,  dirt  and  grass,  coniferous  trees,  and  deciduous  trees  for  the  Gray¬ 
ling  tests;  and  grass,  brittle  bush,  creosote  bush,  mixed  bushes,  and 
saguaro  cactus  for  the  Yuma  tests.  Since  the  trees  are  not,  in  general,  large 
enough  to  fill  an  entire  radar  cell,  the  cells  identified  as  a  coniferous  or  de¬ 
ciduous  tree  may  also  include  some  data  from  grass  or  dirt,  which  will 
have  a  slight  effect  on  the  absolute  RCS  measurements  of  the  various  cell 
types.  We  averaged  the  measured  RCS  in  square  meters  at  each  frequency 
over  all  ramps  of  data  for  each  individual  cell  to  compute  the  cell  average 
for  each  data  run.  We  then  converted  the  RCS  to  the  RCS  per  unit  area  (o^) 
by  dividing  by  the  beam  area  intersecting  the  groimd  at  the  cell  range.  We 
averaged  the  values  and  used  them  to  compute  the  average  of  o®  in  m^/ m^ 
and  the  standard  deviation  in  dB  space  for  each  cell  for  the  entire  measure¬ 
ment  period.  Tables  2  and  3  give  the  composite  RCS  statistics  and  standard 
deviation  for  these  selected  cells  from  the  Grayling  field  tests  for  both  cir¬ 
cular  and  linear  polarization,  respectively. 

The  individual  cell  statistics  shown  in  tables  2  and  3  can  be  used  to  make 
some  general  statements  about  the  cell  RCS's  for  both  circular  and  linear 
polarization  at  the  Grayling  site.  It  is  apparent  that  the  LR  and  RL  RCS  val¬ 
ues  are  about  equal,  and  that  the  RR  and  LL  RCS  values  are  also  about 
equal  for  all  cell  types.  Also,  the  LR  RCS  is  about  4  dB  larger  than  the  RR 


Table  2.  Mean  RCS 

Cell  type 

Location 

RR 

LL 

LR 

RL 

KAmJ  /III  }  dllU 

standard  deviations 

Dirt  and  grass 

Gray  I 

-23.1  + 1.9 

-23.011.9 

-19.211.5 

-19.1  ±  1.5 

in  dB  for  selected 

Grass  and  bushes  Gray  I 

-24.5  ±  2.4 

-24.6 1 2.5 

-20.611.8 

-20.711.8 

cells  for  circular 

Coniferous 

Gray  I 

-22.2  +  1.1 

-22.211.2 

-18.1 1 1.4 

-18.3  +  1.2 

polarization  from 

Deciduous 

Gray  I 

-23.6  ±  3.5 

-23.3 1 3.5 

-16,813.9 

-17.214.0 

Grayling  field  tests. 

Grass  1 

Gray  II 

-24.0  ±  3.8 

-23.514.0 

-21.412.2 

-21.212.2 

Grass  2 

Gray  II 

-24.4  ±  3.3 

-24.3 1 3.4 

-21.211.9 

-21.111.8 

Grass  3 

Gray  II 

-25.3 1  3.6 

-24.6 1 3.7 

-22.9 1 2.1 

-22.412.1 

Grass  4 

Gray  II 

-22.4 1 4.3 

-22.4 1 4.6 

-20.7 1 2.5 

-20.5 1 2.5 

Grass  and 

Gray  II 

-25.1 1  2.0 

-22.2 1 2.2 

-20.6 1 1.4 

-20.411.3 

coniferous 

Coniferous  1 

Gray  II 

-19.0 1 2.0 

-18.712.0 

-15.9 1 1.3 

-15.611.3 

Coniferous  2 

Gray  II 

-22.6 1 2.0 

-20.712.0 

-19.0 1 1.4 

-18.411.3 

Coniferous  3 

Gray  II 

-20.1 1 2.0 

-20.1 1 2.0 

-17.1 1 1.3 

-16.811.3 

Grass  and 

Gray  II 

-22.2 1 3.2 

-21.213.4 

-19.3 1 2.1 

-18.9 1 2.1 

deciduous 

Deciduous  1 

Gray  II 

-20.1 1 2.7 

-19.712.2 

-16.6 1 1.4 

-16.411.3 

Deciduous  2 

Gray  II 

-21.2 1 2.9 

-21.012.6 

-17.611.4 

-17.311.3 

Deciduous  3 

Gray  II 

-23.0 1 2.1 

-22.512.1 

-19.0 1 1.4 

-18.811.4 
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Table  3.  Mean  RCS 
dB  (m^/m^)  and 
standard  deviations 
in  dB  for  selected 
cells  for  linear 
polarization  from 
Grayling  field  tests. 


Cell  type 

Location 

w 

HH 

HV 

VH 

Dirt  and  grass 

Gray  I 

-18.3  ±1.7 

-18.0  ±1.7 

-25.8  ±  2.3 

-25.8  ±  2.4 

Grass  and  bushes 

Gray  I 

-19.3  ±  1.9 

-19.6  ±1.9 

-27.3  ±2.6 

-27.3  ±  2.7 

Coniferous 

Gray  I 

-17.6  ±1.0 

-17.5  ±  1.0 

-24.5  ±1.2 

-24.5  ±1.4 

Deciduous 

Gray  I 

-16.4  ±2.3 

-15.8  ±  2.5 

-25.8  ±  2.2 

-25.9  ±2.4 

Grass  1 

Gray  II 

-20.3  ±  2.5 

-20.2  ±  2.6 

-27.4  ±3.9 

-27.3  ±3.9 

Grass  2 

Gray  II 

-20.1  ±  2.1 

-20.5  ±  2.1 

-28.5  ±  3.2 

-28.5  ±3.3 

Grass  3 

Gray  II 

-21.4  ±2.6 

-21.6  ±2.4 

-29.0  ±3.4 

-28.7  ±3.5 

Grass  4 

Gray  II 

-19.2  ±3.1 

-19.7  ±2.9 

-26.1  ±  4.6 

-26.0  ±  4.6 

Grass  and 

Gray  II 

-21.0  ±1.5 

-18.3  ±  1.3 

-27.4  ±2.0 

-27.4  ±  2.0 

coniferous 

Coniferous  1 

Gray  II 

-15.2  ±1.4 

-15.0  ±1.3 

-22.6  ±1.9 

-22.5  ±  1.8 

Coniferous  2 

Gray  II 

-18.9  ±1.5 

-17.0  ±1.4 

-25.3  ±1.8 

-24.8  ±  1.8 

Coniferous  3 

Gray  II 

-16.5  ±  1.4 

-16.1  ±  1.4 

-23.7  ±2.1 

-23.6  ±  2.0 

Grass  and 

Gray  II 

-18.5  ±  2.4 

-17.9  ±  2.4 

-25.8  ±  3.3 

-25.4  ±  3.3 

deciduous 

Deciduous  1 

Gray  II 

-16.0  ±1.7 

-15.5  ±1.5 

-24.0  ±  2.3 

-23.8  ±  2.3 

Deciduous  2 

Gray  II 

-17.1  ±  1.8 

-16.3  ±1.5 

-25.1  ±  2.7 

-24.9  ±  2.7 

Deciduous  3 

Gray  II 

-18.5  ±  1.6 

-18.0  ±1.5 

-26.7  ±2.0 

-26.5  ±  2.0 

RCS  for  all  the  cells  for  Grayling  I,  except  for  the  deciduous  tree,  where  the 
ratio  (LR/RR)  is  almost  7  dB.  For  Grayling  II,  this  ratio  is  between  2  and  4 
dB,  with  the  average  about  3  dB.  For  linear  polarization,  the  RCS's  for  W 
and  HH  are  about  equal,  and  the  RCS's  for  HV  and  VH  are  almost  identi¬ 
cal.  In  addition,  the  W  RCS  averages  about  7  dB  larger  than  the  HV  RCS 
for  all  cell  types.  The  standard  deviation  values  in  the  tables  indicate  that 
the  RCS  for  the  deciduous  tree  cells  had  the  largest  variation  over  the 
measurement  period  for  Grayling  I;  however,  this  variation  is  due  almost 
entirely  to  a  6-dB  decrease  in  the  average  RCS  of  the  deciduous  tree  cell 
from  the  beginning  of  the  test  to  the  end.  In  general,  the  grass  cells  actually 
had  the  largest  variation  in  RCS  due  to  changing  weather  conditions  for 
both  tests;  for  Grayling  II,  the  deciduous  tree  cells  were  next,  with  the  ever¬ 
green  tree  cells  having  the  smallest  dependency  on  weather. 

To  identify  the  causes  of  the  variation  in  the  measured  RCS  for  individual 
cells,  the  RCS  data  were  plotted  versus  time  and  compared  with  the 
weather  data  observed  at  those  times.  The  behavior  of  the  four  selected  cell 
types  from  Grayling  I  are  shown  in  figure  2,  in  which  the  average  RCS's  for 
each  data  run  for  each  cell  type  are  plotted  versus  run  number.  The  run 
numbers  are  sequential  in  time  from  the  beginning  to  the  end  of  the  meas¬ 
urements.  The  weather  conditions  present  during  the  data  runs  are  shown 
in  figure  3.  Upon  comparison  of  figure  2  with  figure  3,  it  is  apparent  that 
the  major  variations  of  the  grass  cells  are  primarily  due  to  the  presence  of 
water,  frost,  or  snow  on  the  ground.  In  fact,  the  largest  variation  in  the  RCS 
for  these  cells  is  due  to  the  refrozen-snow  condition  observed  for  nm  num¬ 
ber  53.  This  is  shown  clearly  in  figure  4,  where  the  grass  cell  RCS  is  plotted 
on  an  expanded  scale,  showing  the  weather  conditions  responsible  for  the 
significant  variations.  The  RCS's  for  the  coniferous  trees  were  somewhat 
less  sensitive  to  these  conditions,  but  the  effects  were  still  evident.  The  RCS 
for  the  deciduous  tree  seemed  to  be  less  sensitive  to  the  weather,  but  the 
RCS  decreased  as  the  leaves  changed  color  and  lost  moisture.  The  results 
for  linear  polarization  are  similar,  and  will  not  be  discussed  here. 
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Figure  2.  Graphs  of 
sample  cell  RCS  (a®) 
versus  run  number 
for  Grayling  I  test. 
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The  Grayling  II  field-test  data  were  taken  during  the  time  period  from 
early  March  to  mid  April.  During  this  time,  there  was  a  considerable 
amount  of  snow  cover — ^both  fresh  snow  and  refrozen  snow  were  present 
during  the  missions.  The  presence  of  the  snow  had  a  considerable  impact 
on  the  measured  RCS  of  the  grass  cells,  with  slightly  less  of  an  effect  on  the 
deciduous  tree  cells,  and  little  effect  on  the  coniferous  tree  cells.  To  under¬ 
stand  and  compare  the  measurement  results  between  cells  of  the  same 
type  and  between  cells  of  different  types,  the  RCS  data  was  plotted  versus 
mission  number  for  all  the  selected  cells  from  Grayling  II  for  both  circular 
and  linear  polarization.  During  missions  in  which  more  than  one  data  nm 
for  linear  or  circular  polarization  was  measured,  the  data  runs  were  aver¬ 
aged  together  to  get  one  data  point  for  each  polarization  for  each  mission. 
The  RCS  data-versus-mission  number  for  selected  cells  are  shown  in  fig¬ 
ures  5  through  10,  with  the  corresponding  weather  conditions  shown  in 
figure  11.  By  comparing  the  weather  data  with  the  various  RCS  plots  for 
specific  cell  types,  the  following  conclusions  can  be  drawn. 
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Figure  3.  Measured 
snow  depth, 
precipitation  rate,  and 
temperature 
coincident  with 
Grayling  I  data  runs. 
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Figure  4.  The  grass¬ 
cell  RCS  (a^)  plotted 
versus  run  number 
with  the  weather 
conditions 
responsible  for  the 
variations  indicated. 
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Figure  5,  RCS  (a®)  for 
indicated  cells  versus 
mission  number  for 
circular  polarization 
(top)  and  linear 
polarization  (bottom) 
for  Grayling  II  test 
on  grass  1  and  2. 
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Figure  6.  RCS  (a^)  for 
indicated  cells  versus 
mission  number  for 
circular  polarization 
(top)  and  linear 
polarization  (bottom) 
for  Grayling  II  test 
on  grass  3  and  4. 
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Figure  8.  RCS  (a^)  for 
indicated  cells  versus 
mission  number  for 
circular  polarization 
(top)  and  linear 
polarization  (bottom) 
for  Grayling  II  test 
on  coniferous  tree  2 
and  3. 
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Figure  9.  RCS  (a®)  for 
indicated  cells  versus 
mission  number  for 
circular  polarization 
(top)  and  linear 
polarization  (bottom) 
for  Grayling  II  test 
on  grass  and 
deciduous  tree. 
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Figure  10.  RCS  (a®) 
for  indicated  cells 
versus  mission 
number  for  circular 
polarization  (top) 
and  linear 

polarization  (bottom) 
for  Grayling  II  test 
on  deciduous  tree  2 
and  3. 
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Grass  cell  RCS  is  dependent  on  the  amount  of  snow  cover  and  snow  condi¬ 
tion  present.  The  RCS  of  the  grass  cells  was  highest  for  refrozen  snow.  (Re¬ 
frozen  snow  can  be  assumed  to  exist  whenever  there  is  snow  cover  and  the 
temperature  has  risen  above  freezing  and  then  dropped  below  freezing.) 
This  occurred  from  missions  15  to  40  and  at  a  few  other  missions  during 
the  test.  For  fresh  snow,  the  measured  RCS  was  slightly  lower  (missions  92 
to  94,  and  131  to  134).  The  RCS  was  the  lowest  for  melting  snow  when 
large  areas  of  water  accumulated  on  the  ground  (missions  9  to  13),  and 
somewhere  in  between  the  fresh-  and  melting-snow  conditions  for  the 
clear  or  no-snow-cover  condition.  The  polarization  ratio  LR/RR  was  about 
3  dB  for  fresh  snow,  about  0  dB  for  refrozen  snow,  and  about  5  dB  for  no 
snow.  For  linear  polarization,  the  polarization  ratio  VV/HV  was  about 
6  dB  for  fresh  snow,  about  4  dB  for  refrozen  snow,  and  about  8  dB  for  the 
no-snow  condition.  It  was  also  apparent  that,  for  the  condition  without 
snow  cover,  there  is  a  slight  decrease  in  the  grass  RCS  over  time,  due  to  the 
increase  in  the  soil's  moisture  content  that  occurs  as  the  ground  thaws  in 
the  spring.  There  was  no  apparent  growth  in  the  grass  during  the  test,  so 
that  any  effect  due  to  the  grass  growing  was  minimal.  The  graphs  of  the 
linear  and  circular  polarization  exhibit  similar  behavior  over  time,  and  the 
amplitudes  of  the  individual  grass  cells  track  very  well  over  time,  with 
only  slight  differences  in  the  absolute  magnitude  of  the  RCS. 

Upon  examination  of  the  graphs  of  the  RCS  versus  mission  number  for  the 
coniferous  tree  cells,  some  interesting  features  were  evident.  First  of  all,  the 
RCS  amplitude  of  the  individual  cells  tracked  each  other  very  well  over 
time  for  both  linear  and  circular  polarization.  There  does  not  appear  to  be 
any  overall  trend  in  these  measurements,  but  some  slight  effects  due  to 
fresh  snow  are  apparent  (missions  131  to  134),  and  there  are  a  large  num¬ 
ber  of  changes  in  the  RCS  that  cannot  be  explained  by  precipitation  and 
snow  cover.  The  circular  polarization  ratio  generally  runs  about  3  to  4  dB 
for  the  entire  test,  with  a  few  missions  where  it  drops  to  around  0  dB.  The 
linear  polarization  ratio  is  generally  about  7  dB,  but  drops  to  around  3  dB 
when  the  circular  polarization  ratio  drops.  By  correlating  the  coniferous 
tree  RCS  measurements  to  the  temperature,  we  found  that  there  was  a 
measurable  increase  in  the  cell  RCS  for  both  linear  and  circular  polariza¬ 
tion  when  the  temperature  dropped  below  freezing.  This  effect  was  prob¬ 
ably  due  to  the  freezing  of  the  needles  and  the  resulting  change  in  the  in¬ 
dex  of  refraction  caused  by  the  conversion  of  the  liquid  water  in  the 
needles  to  ice.  A  more  in-depth  study  of  this  behavior  will  be  conducted  in 
the  near  future,  but  is  not  included  in  this  report. 

The  graphs  of  the  RCS-versus-mission  number  for  the  deciduous  tree  cells 
exhibit  behavior  similar  to  the  grass  cells.  Since  these  trees  were  not  fully 
foliated  during  the  measurement  period,  the  measured  RCS  is  a  combina¬ 
tion  of  the  ground  around  the  trees  and  the  trees  themselves.  The  magni¬ 
tude  of  this  effect  can  be  seen  by  comparing  figures  5  and  6  to  figures  9 
and  10.  These  figures  track  each  other  quite  well,  especially  during  the 
refrozen-snow  condition,  but  the  total  RCS  for  the  deciduous  trees  is  be¬ 
tween  2  and  3  dB  higher.  This  would  indicate  that  the  ground  is  contribut¬ 
ing  almost  half  of  the  total  RCS.  There  is  even  a  slight  decrease  in  the  RCS 
values  over  time  for  circular  polarization,  as  was  evident  in  the  grass  cells. 


The  polarization  ratios  for  both  linear  and  circular  polarization  are  compa¬ 
rable  to  those  of  the  grass  cells.  Although  it  was  not  observed,  the  RCS  of 
deciduous  trees  should  increase  as  the  leaves  grow  in  a  similar  manner  to 
the  decrease  shown  for  Grayling  I,  when  the  leaves  of  the  deciduous  tree 
turned  brown.  To  separate  the  effects  of  the  ground  return  on  the  total 
RCS,  analysis  is  required  at  a  higher  range  resolution,  so  that  the  contribu¬ 
tion  from  the  ground  can  be  removed  and  the  contribution  from  the  tree 
can  be  isolated.  This  analysis  is  planned  for  the  near  future. 

We  started  the  Yuma  field  test  about  mid-March  and  concluded  it  at  the 
end  of  April  1993.  Due  to  the  excessive  amount  of  rain  that  fell  during  the 
winter,  blooming  of  the  desert  occurred  well  before  the  first  measurement 
mission.  In  addition,  there  was  only  one  rain  event  during  the  entire  test 
period.  In  tables  4  and  5,  the  composite  RCS  statistics  for  selected  cells  over 
the  entire  test  period  are  shown.  The  data  indicate  that  there  was  very  little 
variation  in  the  cell  RCS's  over  the  test  period  for  all  cells  except  the 
saguaro  cactus.  As  with  the  Grayling  data,  the  RL  and  LR  amplitudes  are 
almost  identical,  and  the  RR  and  LL  amplitudes  are  very  similar.  Likewise, 
the  VH  and  HV  returns  are  identical  for  linear  polarization,  and  the  W 
and  HH  returns  are  very  similar.  The  measured  cell  polarization  ratios  are 
also  similar  to  those  measured  in  Grayling  1,  and  slightly  higher  than  those 
observed  in  Grayling  II  for  all  cells  except  the  saguaro  cactus,  where  the 
polarization  ratios  were  6.5  and  16  dB  for  circular  and  linear  polarization, 
respectively.  The  Grayling  II  polarization  ratios  were  slightly  lower  for 
most  cells  due  to  the  effects  of  the  snow  cover  and  variable  weather  condi¬ 
tions,  which  tend  to  increase  the  returns  for  RR  and  LL  compared  to  RL 
and  LR.  To  examine  the  behavior  of  the  cell  RCS  amplitudes  over  time  for 
the  Yuma  test,  the  RCS's  of  the  selected  cells  were  plotted  versus  mission 
number.  These  plots  are  shown  in  figures  12, 13, 14,  and  15  for  both  circu¬ 
lar  and  linear  polarization,  respectively;  corresponding  weather  data  is 
shown  in  figure  16. 


Table  4.  Mean  RCS  (^ell  type  Location  RR  LL  LR  RL 

dB  (m^/m^)  and  - — — - 


standard  deviations 

Brittle  bush  1 

Yuma 

-23.7  ±  0.5 

-23.8  ±  0.6 

-18.9  ±  0.6 

-18.9  ±  0.6 

in  dB  for  selected 

Brittle  bush  2 

Yuma 

-27.6  ±  0.6 

-27.1  ±0.7 

-23.1  ±  0.6 

-23.4  ±  0.6 

cells  for  circular 

Creosote  bush  1 

Yuma 

-25.6  ±  0.7 

-26.1  +  0.7 

-21.8  ±0.9 

-21.9  ±0.9 

polarization  from 

Creosote  bush  2 

Yuma 

-24.1  ±  0.7 

-24.1  ±  0.8 

-19.3  ±1.1 

-19.3  ±1.1 

Yuma  field  test. 

Mixed  bushes  1 

Yuma 

-24.7  ±  0.8 

-26.2  ±  0.8 

-19.3  ±1.0 

-19.4  ±1.0 

Mixed  bushes  2 

Yuma 

-20.1  ±0.8 

-20.0  ±  0.7 

-15.9  ±  0.8 

-15.9  ±  0.8 

Grassy  area 

Yuma 

-29.3  ±  1.9 

-28.5  ±1.4 

-24.9  ±1.1 

-24.6  ±  1.0 

Saguaro  cactus 

Yuma 

-18.1  ±2.1 

-18.0  ±2.1 

-11.6  ±2.1 

-11.6  ±2.1 

Table  5.  Mean  RCS 
dB  (m^/m^)  and 

Cell  type 

Location 

w 

HH 

HV 

VH 

standard  deviations 

Brittle  bush  1 

Yuma 

-17.9  ±  0.4 

-18.1  ±  0.6 

-26.9  ±  0.6 

-26.7  ±0.5 

in  dB  for  selected 

Brittle  bush  2 

Yuma 

-22.2  ±  0.6 

-22.0  ±  0.5 

-30.2  ±  0.5 

-30.3  ±  0.6 

cells  for  linear 

Creosote  bush  1 

Yuma 

-20.5  ±  0.8 

-21.3  ±0.9 

-28.6  ±  0.8 

-28.5  ±  0.8 

polarization  from 

Creosote  bush  2 

Yuma 

-17.9  ±1.1 

-18.2  ±1.0 

-27.6  ±  0.8 

-27.5  ±  0.8 

Yuma  field  test. 

Mixed  bushes  1 

Yuma 

-17.8  ±  0.9 

-19.8  ±  0.8 

-28.2  ±  0.6 

-28.1  ±  0.6 

Mixed  bushes  2 

Yuma 

-14.9  ±  0.6 

-14.9  ±  0.7 

-23.1  ±  0.8 

-23.1  ±  0.8 

Grassy  area 

Yuma 

-24.1  ±  1.3 

-23.7  ±1.0 

-32.6  ±1.6 

-32.1  ±  1.4 

Saguaro  cactus 

Yuma 

-8.7  ±  2.2 

-11.4  ±  2.1 

-27.8  ±2.2 

-27.7  ±  2.1 
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Figure  12.  RCS  (a®) 
for  indicated  cells 
versus  mission 
number  for  circular  o 
polarization  (top)  and  ^ 
linear  polarization  ^ 
(bottom)  for  Yuma 
test  on  brittle  bush  1 
and  2. 
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Figure  13.  RCS  (a^) 
for  indicated  cells 
versus  mission 
number  for  circular 
polarization  (top)  and 
linear  polarization 
(bottom)  for  Yuma 
test  on  creosote  bush 
1  and  2. 
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Figure  14.  RCS  (cfi) 
for  indicated  cells 
versus  mission 
number  for  circular 
polarization  (top)  and 
linear  polarization 
(bottom)  for  Yuma 
test  on  mixed  bushes 
1  and  2. 
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Figure  16.  Measured 
soil-moisture  levels 
averaged  over  the 
ASL  and  TIPS  sites 
and  temperature 
observed  during 
missions  at  Yuma. 
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The  graphs  of  the  cell  RCS's  for  the  bushes  shown  in  figures  12  to  14  show 
very  little  variation  in  the  cell  amplitudes  over  the  measurement  missions. 
The  grassy  area  in  figure  15  shows  a  slight  decrease  in  the  RCS  over  time, 
due  in  part  to  a  decrease  in  the  moisture  content  of  the  grass  as  it  turned 
brown  during  the  test.  The  saguaro  cactus  (also  shown  in  fig.  15)  exhibited 
large  excursions  in  the  measured  RCS  for  both  linear  and  circular 
polarization  over  the  various  missions.  This  effect  was  even  larger  than  in¬ 
dicated  by  the  figure  when  the  individual  data  runs  were  compared  in¬ 
stead  of  the  mission  averages.  At  first,  this  effect  was  puzzling,  but  after 
carefully  examining  the  data,  we  determined  the  cause  of  these  variations: 
the  variations  were  caused  by  the  motion  of  the  cactus  in  the  wind.  When 
there  was  little  or  no  wind,  the  measured  RCS  was  relatively  stable.  How¬ 
ever,  when  the  wind  blew  strongly,  there  were  wild  variations  in  the  meas¬ 
ured  RCS.  The  cactus  is,  essentially,  a  combination  of  cylinders  oriented 
vertically,  and  the  radar  illuminates  the  cactus  almost  perpendicular  to  the 
long  dimension  of  these  cylinders.  These  multiple  cylinders  produce  an  in¬ 
terference  pattern  over  the  bandwidth  of  the  radar  and  this  pattern  varies, 
depending  on  the  absolute  orientation  of  the  cylinders.  This  effect  is  the 
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cause  of  the  very  high  linear  polarization  ratio.  In  addition,  the  cactus  pre¬ 
sents  a  large  aperture  to  the  radar  and,  consequently,  because  of  diffrac¬ 
tion,  generates  a  very  narrow  reflection  pattern.  A  slight  variation  in  the 
angle  of  illumination  of  the  cactus  can  cause  a  large  variation  in  the  inter¬ 
ference  and  the  reflection  pattern  characteristics.  These  effects  occur  when 
the  wind  causes  the  cactus  to  bend  or  to  slightly  change  position  in  the 
wind.  The  weather  data  indicates  a  slight  increase  in  the  soil  moisture  at 
mission  60.  This  change  took  place  as  a  result  of  the  only  rainfall  during 
the  entire  test.  The  soil-moisture  levels  were  very  low  over  the  entire  test 
period,  with  an  overall  decrease  over  time.  By  examining  the  data  care¬ 
fully,  we  found  no  obvious  trends  that  could  be  attributed  to  soil  moisture, 
humidity,  or  temperature.  To  isolate  the  different  environmental  effects,  we 
performed  the  more  in-depth  analysis  described  below. 


3.2  Clutter  Distribution  Analysis 

We  also  used  the  collected  data  to  generate  RCS  histograms  of  sample  clut¬ 
ter  cells  for  RR  and  RL  polarization.  We  compiled  these  histograms  from 
over  1800  individual  RCS  measurements  made  across  the  radar  bandwidth 
for  each  cell  from  the  Grayling  I  and  Yuma  tests,  choosing  the  files  to  give 
a  representative  comparison  of  distributions  for  the  different  clutter  types. 
The  histograms  are  shown  in  figures  17  and  18,  where  the  number  of  meas¬ 
urements  occurring  in  a  1-dB  bin  is  plotted  versus  the  measured  RCS  per 
unit  area  (o'^)  of  the  bin.  It  was  apparent  that  each  type  of  clutter  cell  had 
about  the  same  distribution  function  for  both  polarizations,  with  a  slight 
difference  in  the  mean  value  of  the  distribution.  The  distribution  function 
for  this  type  of  histogram  is  dependent  on  the  number  and  magnitude  of 
the  individual  scatterers  in  the  clutter  cell.  For  instance,  the  distribution  for 
a  single  large  scatterer  (e.g.,  a  comer  reflector)  would  be  narrowly  peaked 
about  the  reflector  RCS  value.  Distributed  clutter  would  have  a  wider  dis¬ 
tribution  similar  to  those  shown  in  the  figures.  Since  the  grass  cells  seem  to 
exhibit  the  most  variation  due  to  changing  weather  conditions,  we  gener¬ 
ated  comparison  histograms  for  the  dirt  and  grass  cells  for  different 
weather  conditions.  These  histograms  are  shown  in  figure  19  for  both  RR 
and  RL  polarizations.  The  figures  show  the  distributions  at  the  beginning 
and  end  of  the  measurements,  for  frost  and  refrozen  snow.  The  histograms 
for  this  particular  cell  were  essentially  the  same  for  the  beginning  and  the 
end  of  the  test,  while  the  frost  condition  modified  the  distribution  only 
slightly.  However,  the  refrozen-snow  condition  distributions  for  both  po¬ 
larizations  became  more  irregular  with  large  peaks,  compared  to  the  be¬ 
ginning  distributions. 

The  distribution  for  both  polarizations  was  almost  the  same  for  this  case, 
as  well.  Analyses  of  the  histograms  for  the  other  cell  types  indicated  that 
the  distributions  did  not  change  significantly  during  the  test.  These  results 
s*Jgg6St  that  a  single  type  of  distribution  function  should  be  able  to  de¬ 
scribe  the  clutter  for  all  conditions,  with  the  exception  of  refrozen  snow. 
Further  analysis  is  needed  to  determine  the  appropriate  functions  for  de¬ 
scribing  the  statistics  of  individual  clutter  cells,  and  for  determining  a  dis¬ 
tribution  function  that  best  describes  the  entire  clutter  area. 
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The  spatial  distribution  of  the  clutter  is  very  important  for  determining 
and  predicting  the  performance  of  any  MMW  system  operating  in  a  clutter 
environment.  To  this  end,  spatial  distributions  have  been  calculated  for 
both  the  Yuma  and  Grayling  II  field  tests  at  both  15-  and  1-m  range  resolu¬ 
tions  for  circular  polarization.  The  spatial  distributions  for  linear  polariza¬ 
tion  are  similar,  and  will  not  be  discussed  here.  The  cross-range  resolution 
of  the  clutter  cells  is  dependent  on  range,  and  varies  from  about  2  to  4  m. 
This  means  that  the  resolution  cells  area  will  be  between  2  and  4  m^  for  the 
high-resolution  case,  and  between  30  to  60  m^  for  the  low-resolution  case. 
These  areas  are  quite  different  in  spatial  extent,  but  both  spatial  ranges  are 
t  important  to  system  operation.  The  low-resolution  spatial  distributions  are 

more  important  for  target  detection,  while  high-resolution  spatial  distribu¬ 
tions  are  more  important  for  tracking.  We  calculated  the  spatial  distribu¬ 
tions,  including  both  the  areas  viewed  by  the  ASL  and  the  TIPS  sensors  at 
Yuma,  and  the  entire  E1/E2  zone  at  Grayling.  We  identified  the  measure¬ 
ment  areas  within  the  radar  scan  pattern  by  overlaying  the  radar  data  on 
the  survey  data.  When  we  made  the  multiple  elevation  measurements  at 
the  same  azimuth  and  range,  we  identified  the  cells  with  the  maximum 
RCS  and  included  them  in  the  final  data  set.  We  used  this  final  set  of  cells, 
which  includes  one  cell  at  each  range  and  azimuth  position  within  the 
measurement  areas,  for  the  spatial  analysis.  A  sample  low-resolution  clut¬ 
ter  map  for  Grayling  II  is  shown  in  figure  20.  For  the  low-resolution  case, 
the  spatial  amplitude  distributions  were  calculated  by  binning  the  RCS 
values  in  1-dB  increments,  and  plotting  the  number  within  each  bin. 
Samples  of  the  low-resolution  spatial  histograms  for  Grayling  and  Yuma 
for  a  clear,  bare-ground  condition  are  shown  in  figures  21  and  22.  These 
distributions  include  all  clutter  types  and  shadowed  cells.  From  the  fig¬ 
ures,  it  is  apparent  that  the  Grayling  clutter  has  a  somewhat  wider  distri¬ 
bution  function  than  the  clutter  at  Yuma,  but  the  overall  shape  and  mean 
of  the  distributions  are  similar.  It  is  also  apparent  that,  for  each  site,  the  RL 
and  LR  distributions  are  similar,  and  the  RR  and  LL  distributions  are  simi¬ 
lar.  The  additional  width  of  the  Grayling  distribution  is  caused  by  the  large 
number  of  cells  shadowed  by  the  trees  at  Grayling  that  do  not  exist  at  the 
Yuma  site. 

To  calculate  the  1-m-resolution  histograms,  we  first  converted  the  data  to 
0.23-m  resolution  by  performing  a  EFT  of  the  frequency  data  for  the  cells 
within  the  regions  that  were  selected  for  spatial  analysis.  For  Yuma,  we  cal¬ 
culated  a  64-pt  FFT,  while  for  Grayling  a  128-pt  FFT  was  needed.  We  con¬ 
verted  each  individual  ramp  to  high  resolution,  and  averaged  the  power  in 
^  the  individual  FFTs  to  give  one  high-resolution  FFT  for  each  cell.  We  added 

four  individual  cells  within  this  FFT  together  to  get  a  1-m-resolution  (pre¬ 
cisely,  0.94-m)  FFT.  We  then  computed  a  1-m-resolution  clutter  map.  (A 
sample  map  for  Grayling  II  is  shown  in  figure  23).  We  binned  these  high- 
range  resolution  images  in  1-dB  steps  to  produce  spatial  distributions  over 
the  same  clutter  area  as  for  the  low-resolution  case.  (These  histograms  are 
shown  in  figures  24  and  25.)  Notice  that  the  RL  and  LR  distributions  are 
similar,  and  that  the  RR  and  LL  distributions  are  also  similar,  as  in  the  case 
of  the  low-resolution  distributions.  However,  a  comparison  of  these  figures 
shows  some  significant  differences.  For  Yuma,  the  distributions  look 
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almost  normal  on  a  logarithmic  scale,  which  means  the  tmderlying  distri¬ 
bution  would  be  log  normal.  However,  the  Grayling  distributions  seem  al¬ 
most  bimodal  for  RR  and  LL  polarization,  and  there  are  increases  in  the 
amplitude  of  the  distributions  for  RL  and  LR  at  the  same  location  as  the 
lower  peak  for  the  RR  and  LL.  As  we  mentioned  earlier,  this  is  caused  by 
the  number  of  cells  that  are  shadowed  by  the  trees.  In  the  high-resolution 
case,  there  are  more  shadowed  cells  by  percentage  than  in  the  low-resolu¬ 
tion  case.  Even  though  the  distributions  for  Grayling  and  Yuma  are  differ¬ 
ent,  the  upper  part  of  the  distributions  are  almost  identical — it  is  the  lower 
sections  that  are  different.  These  distributions  did  not  change  noticeably 
during  the  Yuma  test,  but  for  Grayling  there  were  significant  changes  in 
weather  that  could  and  did  affect  the  distributions. 

To  demonstrate  the  effects  of  weather  changes  on  the  1-m-resolution  distri¬ 
butions,  four  different  conditions  were  selected  for  analysis.  These  condi¬ 
tions  were:  clear  and  dry  (mission  150),  melting  snow  (mission  10),  fresh 
snow  (mission  133),  and  refrozen  snow  (mission  20).  The  1-m-resolution 
spatial  amplitude  distributions  for  these  conditions  for  circular  polariza¬ 
tion  are  shown  in  figure  26.  The  figure  shows  that  spatial  distributions 
changed  noticeably  for  different  ground  conditions.  The  distributions  for 
RR  polarization  are  more  variable  than  those  for  RL  polarization.  The  dis¬ 
tribution  function  for  RR  is  almost  identical  to  the  distribution  for  RL 
polarization  for  both  the  refrozen-  and  fresh-snow  conditions,  the  only  dif¬ 
ference  being  a  slight  shift  in  the  mean.  This  is  not  the  case  for  the  clear  and 
melting-snow  condition,  where  the  distributions  have  the  largest  differ¬ 
ence  between  polarizations.  The  overall  polarization  ratio  RL/RR  for  the 
clutter  distributions  is  the  lowest  for  refrozen  snow,  somewhat  bigger  for 
fresh  snow,  larger  still  for  clear  and  dry,  and  the  largest  for  the  melting- 
snow  condition.  It  is  also  clear  that  one  distribution  function  will  not  de¬ 
scribe  the  1-m  spatial  distributions  for  all  environmental  conditions.  A 
more  detailed  analysis  is  required  to  best  characterize  both  the  15-  and  1-m 
distributions.  This  analysis  will  be  completed  in  the  future. 
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Figure  21.  Sample 
low-resolution  spatial 
histograms  for 
Grayling  II  clutter. 
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Figure  22.  Sample 
low-resolution  spatial 
histograms  for  Yuma 
clutter. 
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Figure  23.  Sample 
high-range  resolution 
RCS  maps  for  RR  and 
RL  polarization  for 
Grayling  11. 
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Figure  24,  Sample 
high-resolution 
spatial  histograms  for 
Grayling  II  clutter. 
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Figure  25.  Sample 
high-resolution 
spatial  histograms  for 
Yuma  clutter. 
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Figure  26.  Sample 
high-resolution 
spatial  histograms  for 
RR  polarization 
above  and  RL 
polarization  below 
for  Grayling  II  for 
different  ground 
conditions. 
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Empirical  Modeling 

4.1  Modeling  Procedures 

Since  we  described  the  radar  in  section  2,  only  specifications  relevant  to  the 
modeling  effort  will  be  repeated  in  this  discussion.  We  limited  the  model¬ 
ing  effort  to  circular  polarization  because  of  the  limited  amount  of  linear 
data  available  for  the  Yuma  test.  The  radar  pulse  width  was  100  ns,  which 
corresponds  to  a  clutter  cell  length  of  15  m.  The  antenna  had  a  two-way 
beamwidth  of  approximately  1.0°,  which  corresponds  to  a  clutter  cell 
width  of  approximately  2.6  m  at  a  range  of  150  m.  The  signal-to-noise  ratio 
was  approximately  20  to  30  dB. 

We  collected  the  ARL  data  sets  analyzed  over  an  approximately  6-  to-7- 
week  period  in  the  late  spring  of  1993  at  Yuma,  and  in  the  early  spring 
of  1994  at  Grayling.  The  Grayling  location  is  similar  to  a  European 
envirorunent,  and  the  Yuma  location  is  a  desert  environment.  We  collected 
data  at  two  rectangular  sites  at  the  Yuma  location,  one  approximately  30  by 
60  m,  and  the  other  about  60  by  60  m,  at  a  range  of  approximately  100  to 
180  m  from  the  radar.  We  collected  data  at  one  site  approximately  100  by 
100  m  at  Grayling,  at  a  range  beginning  at  100  m  from  the  radar,  and 
ending  at  about  200  m  from  the  radar.  The  radar  depression  angles  were 
between  2°  and  5°. 

We  scanned  the  designated  sites  at  two  or  three  different  radar  depression 
angles  (the  angle  incremented  by  1°  after  each  azimuth  scan),  selecting  the 
scan  patterns  to  maximize  the  amount  of  ground-clutter  data  collected  at 
each  site,  while  minimizing  the  acquisition  time.  We  analyzed  the  data 
with  frequency-averaged  o®  values.  We  selected  clutter  cells  for  analysis  by 
measuring  the  radar  at  each  range  and  azimuth  position  that  yielded  the 
highest  average  among  the  different  elevation  angles,  assigning  each 
clutter  cell  a  specific  vegetation  type  by  examining  an  overlay  of  the  veg¬ 
etation  survey  data  provided  by  Waterways  Experiment  Station  (WES)  on 
high-range-resolution  clutter  maps.  A  breakdown  for  the  cell  types,  num¬ 
ber  of  cells  of  each  type,  and  descriptions  of  the  vegetation  types  deter¬ 
mined  for  the  ground  clutter  cells  at  Yuma  and  Grayling  is  given  in  tables 
6  and  7. 

For  the  Grayling  site,  the  tree  vegetation  types  were  described  as  "large"  or 
"small."  This  descriptor  was  based  upon  the  relative  RCS  contribution  of 
the  trees  in  a  clutter  cell,  compared  to  the  RCS  of  the  surrounding  grass 
and  dirt  and  not  upon  the  absolute  magnitude  of  A  difference  of  ap¬ 
proximately  10  dB  or  greater  was  considered  large.  More  coniferous  trees 
were  classified  as  large,  compared  to  deciduous  trees,  because  the  decidu¬ 
ous  trees  had  no  leaves  during  the  test  and,  therefore,  usually  had  a 
smaller  RCS.  We  did  not  analyze  small  coniferous  vegetation  or  small 
mixed  vegetation  because  of  the  few  clutter  cells  with  those  classifications. 

We  divided  the  tests  into  approximately  three  missions  per  day  (1-hour 
measurement  intervals  selected  at  random)  during  the  data-collection  peri¬ 
ods;  within  each  mission  we  made  several  radar  scans  of  the  entire  desig¬ 
nated  areas.  For  the  Yuma  data  set,  we  usually  made  two  or  three  data  runs 


Table  6.  Description 
of  Yuma,  AZ, 
vegetation  types. 


Table  7.  Description 
of  Grayling,  MI, 
vegetation  t)T>es. 


Vegetation 

Number  of 

type 

Description 

cells 

Trees 

palo  verde  trees,  large  creosote  bushes. 

48 

and/or  cactus 

Mixed 

trees  and  bushes  with  grass,  dirt, 
and  gravel 

12 

Shadowed 

low-lying  or 

79 

shadowed  vegetation 

Vegetation 

Number  of 

type 

Description 

cells 

Grass 

grass  and  dirt  with  an  occasional 

103 

small  bush 

Large  deciduous 

large  deciduous  trees 

10 

Small  deciduous 

small  deciduous  trees  with 

20 

grass  and  dirt 

Large  coniferous 

large  coniferous  trees 

22 

Small  coniferous 

small  coniferous  trees  with 

1 

grass  and  dirt 

Large  mixed 

large  and  small  coniferous  and 

8 

deciduous  trees 

Small  mixed 

small  coniferous  and  deciduous  trees 

2 

with  grass,  bushes,  and  dirt 

Shadowed 

low-lying  or  shadowed  vegetation 

126 

with  circular  polarization  for  each  mission.  These  intra-mission  data  were 
determined  to  be  highly  correlated  and,  therefore,  were  averaged.  For  the 
Yuma  test,  we  used  data  from  a  total  of  80  missions  for  the  analysis.  For  the 
Grayling  data  set,  we  usually  made  two  data  runs  with  circular  polariza¬ 
tion  for  each  mission.  Because  of  the  changing  environmental  conditions 
during  the  missions,  the  intra-mission  data  were  not  highly  correlated  and, 
therefore,  we  did  not  average  them.  For  the  Grayling  test,  we  used  data  for 
a  total  of  253  runs  for  analysis. 


4.2  Analysis 

We  developed  a  multivariate  linear  regressive  model  that  described  the 
long-term  trends  in  o®  for  different  vegetation  types  and  polarizations.  We 
analyzed  the  data  by  averaging  the  values  together  (in  dBm^/ m^)  for 
each  vegetation  type  and  polarization  for  each  mission  for  the  Yuma  data 
set,  and  for  each  run  within  a  mission  for  the  Grayling  data  set.  We  exam¬ 
ined  the  effects  of  all  the  relevant  measured  direct  and  indirect  environ¬ 
mental  parameters  on  the  average  value  of  o®  for  each  vegetation  type  and 
polarization.  We  considered  the  following  parameters  to  be  independent 
variables  for  the  Yuma  a®  model:  chronological  time,  time  of  day  relative  to 
sunrise,  relative  humidity,  absolute  humidity,  temperature,  solar  flux, 
wind  speed,  and  average  soil  moisture  content.  We  measured  relative 
humidity,  temperature,  solar  flux,  and  wind  speed  every  hour,  and  soil 
moisture  content  once  per  day  at  several  locations.  For  the  Grayling  o® 


f 
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model,  we  considered  the  following  additional  parameters:  snow  cover, 
water  cover,  snow  depth,  snow  density,  snow  grain  diameter,  and  the 
equivalent  water  content  of  the  snow.  We  usually  measured  these  param¬ 
eters  once  per  day,  and  also  examined  interactions  between  selected  pa¬ 
rameters.  We  did  not  include  chronological  time  in  the  Grayling  o®  model, 
because  we  surmised  that  it  was  an  indirect  variable  that  could  be  ex¬ 
pressed  through  other  more  direct  variables.  We  examined  all  the  param¬ 
eters  described  above  as  possible  independent  variables  for  the  model  by 
visual  inspections  of  the  plots  of  the  averaged  cfi  values  versus  selected 
parameters,  and  by  analyzing  the  F  values  and  correlation  coefficients  ob¬ 
tained  from  various  least-square-fit  analyses  of  the  data. 

Initially,  we  performed  the  model  calculations  separately  on  o®  values 
measured  with  RL,  LR,  RR,  and  LL  polarization.  The  results  indicated  that 
the  cross-polarization  (RL  and  LR)  and  the  co-polarization  (RR  and  LL)  re¬ 
sults  were  within  the  expected  error  rate.  Although  a  paired-t  test  indi¬ 
cated  that  neither  the  co-  nor  cross-polarization  measurements  were  equal, 
the  differences  were  usually  very  small  (less  than  0.3  dB).  This  result  justi¬ 
fied  averaging  together  the  co-polarization  a®  measured  values  and  aver¬ 
aging  together  the  cross-polarization  o®  measured  values  before  calculat¬ 
ing  the  model. 

4.2.1  Yuma  Data 

Our  analysis  of  the  data  collected  at  Yuma  indicated  that  chronological 
time  was  the  most  significant  parameter,  and  that  soil  moisture  content 
and  relative  humidity  were  the  next  most  significant  parameters.  The  aver¬ 
age  a®  values  measured  for  different  vegetation  types  and  polarizations  at 
Yuma  were  modeled  by 

^ij  -  Poi;  +  +  Psi/^m  +  (1) 

where  ct®  is  in  dBm^/m^,  i  =  vegetation  type,  j  =  co-  or  cross-circular- 
polarization,  P  =  model  coefficients,  T  =  chronological  time  in  hours/ 
10000,  M  =  soil  moisture  content  in  percent/100,  R  =  relative  humidity  in 
percent/ 100,  m  =  index  of  the  corresponding  variable  associated  with  a 
particular  mission  number,  and  e  =  randomly  distributed  error  with  zero 
mean  and  a  standard  deviation  in  dB.  Figures  27  and  28  show  a  plot  of  the 
measured  o®  values  and  the  fit  resulting  from  applying  equation  (1)  with 
the  coefficients  listed  in  table  8.  The  symbols  show  the  measured  data 
points  and  the  lines  show  the  model  calculations.  The  residue  of  the 
modeled  minus  the  measured  a®  values  was  approximately  normally 
distributed. 

The  model  coefficients  and  their  standard  error  rates  for  different  polariza¬ 
tion  and  vegetation  types  are  shown  in  table  8.  You  can  determine  the  sig¬ 
nificance  of  the  coefficients  by  multiplying  the  coefficients  by  the  range  of 
the  corresponding  independent  variables.  The  chronological  time  varied 
from  the  start  of  the  measurements  on  15  March  to  the  end  of  the  measure¬ 
ments  on  30  April,  or  from  0  to  1100  hours.  Soil  moisture  content  varied 
from  0.96  to  5.18  percent,  with  a  mean  of  1.97  percent;  relative  humidity 
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Figure  27.  Model 
calculations  and 
measurements  of  o® 
for  tree  and  grass 
vegetation  at  Yuma, 
AZ,  with  circular 
cross-polarization. 


Figure  28.  Model 
calculations  and 
measurements  of  o® 
for  tree  and  grass 
vegetation  at  Yuma, 
AZ,  with  circular 
co-polarization. 


Table  8.  Linear 
model  coefficients 
calculated  from 
Yuma  o®  values. 


Vegetation 

type 


Polarization  offset 


Chronological 

time 


Soil 

moisture 


Relative 

humidity 


Trees 

Grass 

Mixed 

Trees 

Grass 

Mixed 


cross 

cross 

cross 

CO 

CO 

CO 


-20.1610.09 
-21.6710.12 
-21.8110.14 
-24.25 1 0.06 
-24.69 1 0.08 
-24.85 1 0.10 


3.010.7 
3.811.0 
6.711.1 
-3.1 1 0.5 
-13. 1 0.7 
-9.8 1 0.8 


-1113 

1414 

-1814 

-512 

-1113 

-1013 


-0.510.2 
-0.510.2 
-0.4 1 0.2 
-0.8 1 0.1 
-0.910.1 
-1.010.2 


varied  front  2  to  74  percent,  with  a  mean  of  26  percent.  The  most  significant 
change  in  was  a  1.5-dB  decrease  over  chronological  time  for  the  a®  of 
grass  measured  with  co-polarization.  The  effects  of  soil  moisture  content 
and  relative  humidity  on  varied,  but  on  average  they  were  responsible 
for  a  maximum  change  of  about  0.5  dB.  The  mixed  vegetation  coefficients 
showed  some  unexpected  variations.  Because  mixed  vegetation  is  a  combi¬ 
nation  of  tree  and  grass  vegetation,  we  expected  that  the  mixed  model  co¬ 
efficients  would  fall  between  the  grass  and  tree  coefficients;  however,  this 
did  not  occur  for  all  coefficients,  due  perhaps  to  the  small  sample  size  of 
the  mixed  vegetation  type,  or  the  fact  that  only  small  trees  were  included 
in  this  vegetation  type  (small  trees  may  have  different  clutter  characteris¬ 
tics  than  large  trees). 
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The  most  significant  environmental  parameter  determined  by  the  model 
was  chronological  time,  which  physically  appears  to  describe  the  accumu¬ 
lated  drying  effect  of  the  sun,  primarily  on  grass  vegetation  over  time.  This 
explanation  is  supported  by  the  previous  weather  conditions.  Heavy  rains 
occurred  earlier  than  expected  at  Yuma,  which  resulted  in  the  peak  of  the 
desert  bloom  occurring  before  the  beginning  of  the  test.  The  vegetation  ap¬ 
peared  to  have  had  limited  growth,  and  the  grass  appeared  to  dry 
throughout  the  test.  A  heavy  rain  occurred  between  missions  25  and  26, 
but  because  of  scheduling  problems,  no  measurements  were  made  until 
two  days  later.  The  effect  of  the  rain  was  still  significant,  but  the  changes  in 
o®  probably  would  have  been  larger  if  we  had  made  measurements  closer 
to  the  rain  event.  This  factor  could  have  resulted  in  the  model  attributing  a 
larger  change  in  o®  to  the  soil  moisture  content  «md  relative  humidity  pa¬ 
rameters.  Only  the  cross-polarization  coefficient  for  grass  indicated  an  in¬ 
crease  in  for  an  increase  in  soil  moisture  content.  A  possible  physical  ex¬ 
planation  for  this  observation  is  that  the  grass  absorbed  the  moisture  faster 
than  the  other  vegetation  types.  This  event  could  change  the  transmission 
and  reflection  coefficients  of  the  grass  (rather  than  simply  attenuate  the  ra¬ 
diation),  thereby  increasing  single-bounce,  cross-polarization  scattering, 
while  reducing  multiple-bounce,  co-polarization  scattering. 

The  effectiveness  of  the  model  can  be  evaluated  by  examining  the  results 
in  table  9.  Column  three  shows  the  estimated  standard  deviation  of  the 
modeled  minus  measured  values  in  dB,  and  column  four  shows  the  cor¬ 
relation  coefficient  adjusted  for  the  number  of  degrees  of  freedom  calcu¬ 
lated  for  measured  and  modeled  values.  The  standard  deviation  deter¬ 
mines  a  quantitative  error  between  modeled  and  measured  values,  and 
the  correlation  coefficient  determines  how  well  the  modeled  and  measured 
a®  values  tracked.  A  low  standard  deviation  and  a  high  correlation  indicate 
a  good  fit. 

Table  9  shows  that  the  correlation  between  the  measured  and  modeled  o® 
values  ranged  from  a  high  of  0.91  for  co-polarization  grass  to  a  low  of  0.72 
for  co-polarization  trees.  These  values  indicate  that  the  major  environmen¬ 
tal  parameters  were  identified  and  successfully  used  by  the  model  to  de¬ 
scribe  the  variation  in  cfi  over  the  duration  of  the  measurement  period.  The 
standard  deviation  varied  between  0.13  and  0.29  dB,  which  is  slightly 
larger  than  the  expected  radar  calibration  error  in  a®  of  0.1  dB.  This  indi¬ 
cates  that  another  variable  may  be  required  in  the  model,  or  that  the  veg¬ 
etation  types  may  need  to  be  more  accurately  described.  Because  of  the 
small  number  of  errors  in  the  model  calculations  and  the  limited  number 
of  mission  sample  points,  we  found  that  further  analysis  of  the  data  to  de¬ 
termine  other  model  variables  was  unnecessary. 


Table  9.  Agreement 
between  modeled 
and  measured  a® 

Vegetation 

type 

Polarization 

Standard 
deviation  (dB) 

Correlation 

coefficient 

values  for  Yuma 

Trees 

cross 

0.19 

0.76 

data  set. 

Grass 

cross 

0.25 

0.73 

Mixed 

cross 

0.29 

0.77 

Trees 

CO 

0.12 

0.72 

Grass 

CO 

0.17 

0.91 

Mixed 

CO 

0.21 

0.80 

49 


4.2.2  Grayling  Data 

We  selected  parameters  for  modeling  the  Grayling  data  set  with  the  same 
methodology  that  was  used  for  Yuma;  however,  for  Grayling  we  applied 
the  procedure  separately  to  data  collected  over  similar  environmental  con¬ 
ditions,  rather  than  to  the  entire  data  set.  We  analyzed  five  environmental 
conditions:  refrozen  ground  and  snow,  drying  grotmd,  melting  snow,  tran¬ 
sitional  snow,  and  fresh  snow.  We  identified,  but  did  not  test,  Ihe  following 
environmental  conditions  because  of  lack  of  sufficient  repeat  conditions: 
rain,  freezing  rain,  snow,  wet  snow,  and  water-saturated  ground.  We  as¬ 
signed  each  data  nm  an  environmental  condition  determined  by  a  qualita¬ 
tive  evaluation  of  the  past  and  present  environmental  data.  The  analysis 
indicated  that  refrozen  ground  and  snow,  as  well  as  drying  groimd  could 
be  empirically  modeled,  and  that  melting,  transitional,  and  fresh  snow 
could  not  be  empirically  modeled  with  the  available  environmental  data. 

The  refrozen  ground  and  snow  environmental  condition  consisted  of  refro¬ 
zen  snow  (which  was  often  present  during  the  beginning  of  the  test),  fro¬ 
zen  or  frost-covered  ground  (which  was  often  present  during  the  end  of 
the  test),  and  a  mixture  of  ground  and  snow  conditions  (which  often  oc¬ 
curred  during  the  middle  of  the  test).  The  air  temperature  was  usually  be¬ 
low  0°C.  The  drying-ground  condition  occurred  primarily  during  the 
middle  and  end  of  the  test;  the  air  temperature  was  usually  above  0°C. 
Drying  ground  could,  potentially,  have  some  snow  cover,  but  not  a  large 
amoxmt.  A  soil  moisture  measurement  was  usually  performed  once  a  day 
for  drying  ground. 

Our  analysis  of  the  Grayling  data  set  during  conditions  of  refrozen  groimd 
and  snow  indicated  that  the  most  significant  environmental  parameters 
were  snow  cover  and  air  temperature.  For  drying-ground  conditions,  the 
most  significant  parameters  were  snow  cover,  soil  moisture  content,  abso¬ 
lute  humidity,  and  wind  speed.  We  measured  these  variables  directly  (ex¬ 
cept  for  snow  cover,  which  we  determined  by  examining  snow  maps  gen¬ 
erated  by  WES).  We  computed  the  snow  cover  ratio  by  calculating  the  ratio 
of  the  area  covered  by  snow  to  the  total  area  of  the  site;  this  calculation  was 
often  interpolated  between  measurements.  We  evaluated  the  model  with 
several  permutations  of  the  snow-cover  parameter:  snow-cover  ratio,  log 
of  the  snow-cover  ratio,  and  square  root  of  the  snow-cover  ratio.  The 
square  root  of  the  snow-cover  ratio  parameter  consistently  had  the  best  re¬ 
sults.  The  a®  values  measured  during  the  refrozen-ground  and  snow  con- 
dirion  were  modeled  by 

+  +  +  (2) 

where  is  in  dBm^/ m^,  i  =  vegetation  type,  j  =  co-  or  cross-circular- 
polarization,  P  =  model  coefficients,  S  =  square  root  of  the  snow  cover-ratio 
in  (m^/m^),  t  =  temperature  in  °C,  r  =  index  of  corresponding  variable  as¬ 
sociated  with  a  particular  data  run  number,  and  e  =  randomly  distributed 
error  with  zero  mean  and  a  standard  deviation  in  dB.  The  values  meas¬ 
ured  during  the  drying-ground  condition  were  modeled  by 

=  Poi;  +  Pli/^r  +  Pay-^r  +  Psi/^r  +  P4iy^r  +  ^ijr  (3) 
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Figure  29.  Model 
calculations  and 
measurements  of  cfi 
for  grass  at  Grayling, 
MI,  with  circular 
cross-polarization  for 
the  refrozen-ground 
and  snow  condition. 


Figure  30.  Model 
calculations  and 
measurements  of  a® 
for  grass  at  Grayling, 
Ml,  with  circular  co¬ 
polarization  for  the 
refrozen-ground  and 
snow  condition. 


Figure  31.  Model 
calculations  and 
measurements  of  a® 
for  coniferous  trees  at 
Grayling,  Ml,  with 
circular  cross¬ 
polarization  for  the 
refrozen-ground  and 
snow  condition. 


where  A  =  absolute  humidity  in  g/m^,  M  =  soil  moisture  content  in 
percent/100,  and  W  =  wind  speed  in  m/ s.  Figures  29  to  34  show  the  aver¬ 
aged  values  and  the  fit  resulting  from  applying  equations  (2)  and  (3) 
with  the  coefficients  listed  in  tables  10  and  11. 


51 


Figure  32.  Model  -18 

calculations  and 
measurements  of  a®  -20 

for  coniferous  trees  ^ 

at  Grayling,  MI, 
with  circular  m 

^  —24 

co-polarization  for  ^ 

the  refrozen-ground 
and  snow  condition. 

-28 
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Figure  33.  Model 
calcuations  and 
measurements  of  a® 
for  coniferous  trees  at 
Grayling,  MI,  with 
circular  cross¬ 
polarization  for  the 
drying-ground 
condition. 


Figure  34.  Model 
calcuations  and 
measurements  of 
for  coniferous  trees 
at  Grayling,  Ml, 
with  circular 
co-polarization  for 
the  drying-ground 
condition. 
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The  model  coefficients  and  their  standard  errors  for  refrozen  ground  and 
snow  and  for  drying  ground  are  shown  in  tables  10  and  11.  The  signifi¬ 
cance  of  the  linear  coefficients  on  o®  can  be  determined  by  multiplying  the 
coefficients  by  the  range  of  the  corresponding  independent  variables.  We 
analyzed  data  runs  from  4  March  to  14  April.  During  the  refrozen-ground 
condition,  the  temperature  varied  from  -18.4  to  0.2°C,  with  an  average 
temperature  of-5.5°C.  The  square  root  of  the  snow-cover  ratio  varied  from 
0.0  to  1.00,  with  an  average  of  0.45.  During  the  drying-ground  condition, 
the  soil  moisture  varied  from  4.6  to  43.3  percent  (averaging  22.6  percent); 
absolute  humidity  varied  from  1.7  to  6.1  g/m^  (averaging  3.4  g/m^);  the 
square  root  of  the  snow-cover  ratio  varied  from  0.0  to  0.45  (averaging  0.13); 
and  the  wind  speed  varied  from  1.0  to  5.4  m/s  (averaging  2.6  m/s). 

The  results  for  the  model  during  the  refrozen-groimd  and  snow  condition 
indicate  that,  for  co-polarization,  the  o®  value  of  coniferous  trees  was 
mainly  dependent  on  the  temperature;  the  value  of  grass  was  mainly 
dependent  on  the  snow  cover;  and  the  o®  values  for  the  remaining  vegeta¬ 
tion  types  were  dependent  on  a  combination  of  both  parameters.  The 
model  calculation  of  for  coniferous  trees  for  co-polarization  changed  by 
5.8  dB  over  the  temperature  range,  and  the  model  calculation  of  o®  for 
grass  for  co-polarization  changed  by  8.8  dB  over  the  snow-cover  range. 
The  cross-polarization  calculations  showed  a  similar  trend,  but  with 
smaller  changes  in  a*^.  The  results  for  the  drying-ground  condition  indicate 
that  the  co-  and  cross-polarization  results  for  all  vegetation  types  were 


Table  10.  Linear 
model  coefficients 

Vegetation  type 

Polarization 

offset 

Temperature 

Snow  cover 

calculated  for  the 

Grass 

cross 

-22.44  ±0.14 

-0.0210.02 

4.010.2 

refrozen-ground  and 

Deciduous  large 

cross 

-21.5610.10 

-0.03 1 0.01 

1.510.1 

snow  condition. 

Deciduous  small 

cross 

-22.35  1  0.11 

-0.0110.01 

2.910.2 

Coniferous  large 

cross 

-20.81 1 0.09 

-0.1810.01 

0.9 1 0.1 

Mixed  large 

cross 

-18.98 1 0.09 

-0.09 1 0.01 

1.2  ±0.1 

Grass 

CO 

-27.55  1  0.28 

-0.01 1 0.03 

8.810.4 

Deciduous  large 

co 

-27.77  1  0.16 

-0.0710.02 

4.210.2 

Deciduous  small 

CO 

-27.76  1  0.24 

-0.0210.03 

7.210.4 

Coniferous  large 

CO 

-24.42  1  0.14 

-0.31 1 0.02 

2.010.2 

Mixed  large 

CO 

-23.67  1  0.13 

-0.1910.02 

2.910.2 

Table  11.  Linear  model  coefficients  calculated  for  the  drying-ground  condition. 


Vegetation 

type 

Polarization 

o®  offset 

Soil 

moisture 

Absolute 

humidity 

Snow 

cover 

Wind 

speed 

Grass 

cross 

-22.20 1 0.14 

-1.010.4 

-0.1710.03 

1.910.2 

0.1410.03 

Deciduous  large 

cross 

-20.8710.18 

-0.410.5 

-0.12 1 0.06 

0.510.3 

0.01 1 0.04 

Deciduous  small 

cross 

-22.11 1 0.15 

-1.5 1 0.4 

-0.1510.03 

1.010.2 

0.1310.03 

Coniferous  large 

cross 

-20.54  1  0.15 

-0.210.4 

-0.0610.02 

-0.510.2 

0.1910.03 

Mixed  large 

cross 

-18.8910.20 

+0.5 1 0.5 

-0.13 1 0.05 

0.410.3 

0.1010.05 

Grass 

CO 

-27.13 1 0.17 

-1.110.4 

-0.2410.04 

2.1 1 0.3 

0.2010.04 

Deciduous  large 

CO 

-27.33  1  0.15 

-0.610.4 

-0.2010.03 

1.010.2 

0.05 1 0.03 

Deciduous  small 

CO 

-27.80  1  0.25 

-1.210.6 

-0.19 1 0.06 

1.510.4 

0.2010.06 

Coniferous  large 

CO 

-24.45 1 0.13 

-0.8 1 0.3 

-0.12 1 0.03 

-0.010.2 

0.2410.03 

Mixed  large 

CO 

-23.7410.24 

-0.5 1 0.6 

-0.23 1 0.06 

0.6  1 0.4 

0.1810.06 
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dependent  on  soil  moisture,  absolute  humidity,  snow  cover,  and  wind 
speed,  with  no  one  dominant  parameter.  These  variables  were  responsible 
for  a  maximum  deviation  in  a®  of  1.1  dB,  and  were  often  responsible  for  a 
maximum  deviation  in  a®  of  0.5  dB  or  less,  values  were  usually  reduced 
when  the  parameters  soil  moisture,  absolute  humidity,  and  wind  speed 
indicated  a  reduction  in  the  amoimt  of  water  in  the  groxmd  clutter. 

The  effectiveness  of  the  model  can  be  evaluated  by  examining  the  correla-  j 

tion  coefficients  and  the  standard  deviation  results  shown  in  tables  12  and 
13.  These  tables  have  the  same  format  as  the  tables  given  in  the  discussion 
of  the  Yuma  results.  The  refrozen-ground  and  snow  condition  results  ^ 

shown  in  table  12  indicate  that  the  modeled  and  measured  o®  values  had 
either  a  high  correlation  or  a  small  standard  deviation,  but  not  both.  The 
lowest  correlation  was  0.75  for  large  deciduous  trees  measured  with  cross¬ 
polarization,  along  with  a  corresponding  good  standard  deviation  of 
0.46  dB.  The  highest  correlation  values  were  over  0.90,  but  usually  had 
standard  deviations  over  1.0  dB.  These  results  indicate  that  another  vari¬ 
able  may  be  required  in  the  model.  The  standard  deviation  results  were 
much  better  for  the  drying-groimd  model  results  in  table  13,  but  the  corre¬ 
lation  values  were  slightly  lower.  Overall,  the  grass  and  coniferous  tree 
vegetation  had  better  model  results  than  deciduous  tree  vegetation. 


Table  12.  Agreement 
between  modeled 
and  measured  a® 

Vegetation 

type 

Polarization 

Standard 
deviation  (dB) 

Correlation 

coefficient 

values  for  refrozen- 

Grass 

cross 

0.65 

0.89 

ground  and  snow 

Deciduous  large 

cross 

0.46 

0.75 

condition. 

Deciduous  small 

cross 

0.54 

0.87 

Coniferous  large 

cross 

0.44 

0.89 

Mixed  large 

cross 

0.42 

0.81 

Grass 

CO 

1.32 

0.91 

Deciduous  large 

CO 

0.75 

0.89 

Deciduous  small 

CO 

1.10 

0.91 

Coniferous  large 

CO 

0.64 

0.92 

Mixed  large 

CO 

0.60 

0.91 

Table  13.  Agreement 

between  modeled  and  Vegetation 

measured  values 

Polarization 

Standard 
deviation  (dB) 

Correlation 

coefficient 

for  the  drying-ground  Grass 

cross 

0.23 

0.87 

condition.  Deciduous  large 

cross 

0.29 

0.45 

Deciduous  small 

cross 

0.24 

0.80 

Coniferous  large 

cross 

0.23 

0.64 

Mixed  large 

cross 

0.32 

0.44 

Grass 

CO 

0.27 

0.88 

Deciduous  large 

CO 

0.24 

0.75 

Deciduous  small 

CO 

0.40 

0.71 

Coniferous  large 

CO 

0.21 

0.81 

Mixed  large 

CO 

0.39 

0.65 
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The  radar  calibration  error  varied  for  different  weather  conditions.  For 
warm  weather  conditions,  the  average  calibration  error  for  o®  was  ap¬ 
proximately  0.1  dB,  and  for  snowing  and  freezing  weather  conditions,  the 
calibration  error  was  approximately  0.3  dB.  Snow  and  ice  collected  in  some 
of  the  calibration  reflectors  for  several  runs  during  the  refrozen-groxmd 
and  snow  environmental  conditions.  (We  took  steps  to  minimize  this  ef¬ 
fect,  but  it  could  not  be  entirely  eliminated.)  The  standard  deviation  results 
for  several  vegetation  types  and  polarizations  were  significantly  larger 
than  the  radar  calibration  error. 

The  refrozen-ground  and  snow  model  results,  in  general,  had  poor  quanti¬ 
tative  agreement  with  the  measured  data  relative  to  the  expected 
calibration  error.  The  high-correlation  coefficients  indicated  that  the  major 
parameters  were  included  in  the  model,  but  some  minor  parameters  were 
not  identified.  We  did  not  include  snow  liquid  water  content,  snow  den¬ 
sity,  and  grain-size  diameter  in  the  model,  because  at  the  once-per-day 
sample  rate,  they  were  not  statistically  significant  parameters.  However,  if 
these  parameters  were  sampled  more  often  or  easily  interpolated  between 
measurements,  we  probably  would  have  included  them  in  the  refrozen- 
ground  and  snow  model,  as  well  as  in  the  models  for  the  other  environ¬ 
mental  conditions  that  were  unsuccessfully  modeled. 

The  drying-ground  model  results  had  good  quantitative  agreement  with 
the  measured  cy®  values.  The  standard  deviation  results  ranged  from  0.2  to 
0.4  dB  for  both  polarizations,  but  the  correlation  coefficients  were  low  for 
the  large  deciduous  and  mixed  trees  for  the  cross-polarization  measure¬ 
ments.  A  possible  explanation  for  this  event  is  that  these  vegetation  types 
were  not  affected  as  much  by  the  environmental  conditions,  which  focused 
more  on  characterizing  the  ground.  We  characterized  the  remaining  veg¬ 
etation  types  reasonably  well  by  the  model  for  both  polarizations. 
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5.  Conclusions 


In  this  paper,  we  described  in  detail  the  95-GHz  clutter  measurements 
taken  on  extensively  instrumented  clutter  patches  at  Grayling,  MI,  and 
Yuma,  AZ.  We  made  the  clutter  measurements  during  the  summer-to-fall 
and  winter-to-spring  transitions  at  Grayling,  and  the  early-to-late-spring 
period  at  Yuma.  The  Grayling  areas  that  we  measured  are  similar  to  a  Eu¬ 
ropean  environment,  and  the  Yuma  test  site  is  in  a  desert  environment. 
During  the  measurement  periods,  we  measured  the  same  clutter  area  re¬ 
peatedly,  using  a  random-sampling  plan,  for  both  circular  and  linear  polar¬ 
ization.  The  resulting  database  is  useful  for  evaluating  clutter  statistics  for 
the  two  environments,  and  for  making  the  statistical  comparisons  neces¬ 
sary  for  evaluating  the  SWOE  MMW  clutter  model. 

We  presented  the  RCS  statistics  for  some  selected  clutter  cells  covering  the 
entire  measurement  period  at  each  field  test.  The  cross-polarized  return 
RCS's  RL  and  LR  for  circular  and  VH  and  HV  for  linear  polarization  meas¬ 
ured  within  individual  clutters  cells  are  virtually  the  same,  while  the  co¬ 
polarized  returns  for  both  circular  and  linear  polarization  are  very  similar. 
The  results  for  the  cross-polarized  radar  returns  are,  therefore,  in  excellent 
agreement  with  scattering  theory. 

We  showed  that  the  measured  RCS's  varied  greatly  for  different  environ¬ 
mental  conditions  for  the  Grayling  field  tests.  The  grass-clutter  cells  were 
the  most  sensitive  to  changing  environmental  conditions  (e.g.,  frost,  fresh 
snow,  refrozen  snow,  and  melting  snow),  while  the  RCS  changes  for  conif¬ 
erous  trees  due  to  these  conditions  were  not  as  significant.  The  RCS's  for 
coniferous  trees  were  found  to  increase  noticeably  when  the  temperature 
dropped  below  freezing.  We  also  learned  that  the  deciduous  tree  cells  were 
less  sensitive  to  changes  in  the  weather  while  fully  foliated  than  to  changes 
in  leaf  condition.  During  the  spring  when  leaves  were  not  present,  the  RCS 
of  the  deciduous  trees  fluctuated  almost  as  much  as  the  RCS  for  the  grass 
cells,  due  to  changes  in  the  environmental  conditions.  We  foimd  that  clut¬ 
ter  cells  of  the  same  type  tracked  each  other  very  well.  We  discovered  that 
the  measured  RCS's  for  Yuma  clutter  types  varied  very  little  over  the 
measurement  period,  due  partly  to  the  lack  of  significant  weather  events 
during  that  timeframe. 

We  presented  three  different  types  of  distribution  functions  for  selected 
clutter  cells  for  each  clutter  site.  We  computed  histograms  within  indi¬ 
vidual  cells  across  the  radar  bandwidth,  and  presented  spatial  amplitude 
histograms  at  both  15-  and  1-m  range  resolutions.  We  showed  that  the  fre¬ 
quency  histograms  were  similar  for  all  clutter  types,  with  only  slight  dif¬ 
ferences  in  the  means.  (These  distributions  were  similar  to  a  log-normal 
distribution.)  The  spatial  histograms  were  slightly  different  between  the 
two  clutter  areas:  the  Grayling  15-m  spatial  histograms  were  slightly 
broader  than  those  from  Yuma.  We  also  found  that  the  Grayling  1-m  spa¬ 
tial  distributions  were  almost  bimodal,  whereas  the  Yuma  1-m  spatial  dis¬ 
tributions  indicated  a  log-normal  distribution.  These  differences  were  due 
mainly  to  the  larger  number  of  shadowed  cells  produced  by  the  trees 
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present  in  the  Grayling  clutter  area.  We  found  the  distributions  to  be 
relatively  stable  during  the  Yuma  test,  while  for  Grayling  there  were  large 
variations  caused  by  weather  changes.  The  frequency  and  spatial  distribu¬ 
tions  for  RR  and  LL  polarization  were  more  sensitive  to  the  varying  envi¬ 
ronmental  conditions  than  the  RL  and  LR  distributions. 

The  0®  values  of  different  vegetation  types  measured  with  circular  polar¬ 
ization  at  Yuma  in  the  spring  of  1993  did  not  change  very  much.  The  most 
significant  environmental  parameter,  as  determined  by  an  empirically  cal¬ 
culated  linear  regressive  model  based  upon  least-squared  estimation  of  0°, 
was  chronological  time.  The  physical  significance  of  this  parameter  was 
interpreted  as  the  accumulated  drying  effects  of  the  sim.  The  model  deter¬ 
mined  that  soil-moisture  content  and  relative  humidity  also  had  a  small  ef¬ 
fect  on  0®.  These  parameters  were  significant  in  describing  the  effects 
caused  by  one  large  rain  event. 

The  0®  values  of  different  vegetation  types  measured  with  circular  polar¬ 
ization  at  Grayling  in  the  spring  of  1994  were  highly  dependent  on  envi¬ 
ronmental  conditions.  A  linear  regressive  model  successfully  described  0® 
measurements  made  during  the  environmental  conditions  of  refrozen 
ground  and  snow,  as  well  as  drying  ground;  however,  it  was  unsuccessful 
in  describing  0®  measurements  made  during  the  environmental  conditions 
of  melting,  transitional,  and  fresh  snow.  The  negative  results  were  prob¬ 
ably  caused  by  the  model's  attempt  to  use  important  environmental  pa¬ 
rameters  (e.g.,  snow  density,  grain-size  diameter,  and  liquid-water  content 
of  the  snow)  that  were  sampled  at  a  rate  of  once  a  day.  This  sample  rate 
was  too  low  for  the  model,  and  a  linear  interpolation  between  measure¬ 
ments  was  not  reasonable.  The  0*^  values  measured  during  the  refrozen- 
ground  and  snow  conditions  were  dependent  on  the  parameters  of  snow 
cover  and  air  temperature,  while  0°  for  grass  was  most  dependent  on 
snow  cover,  and  the  0°  value  for  coniferous  trees  was  primarily  dependent 
on  temperature.  The  0*^  values  measured  during  the  drying-ground  condi¬ 
tion  were  dependent  on  the  parameters  of  snow  cover,  soil-moisture  con¬ 
tent,  absolute  humidity,  and  wind  speed,  with  no  one  parameter  having  a 
dominant  effect  for  any  vegetation  t^e.  The  differences  between  the  mea¬ 
sured  and  modeled  0*^  values  were  usually  larger  than  the  expected  cali¬ 
bration  errors  of  the  radar  system,  but  the  correlation  between  the  mod¬ 
eled  and  measured  0®  values  was  usually  high.  This  result  indicates  that, 
for  the  entire  Yuma  data  set  and  for  a  large  portion  of  the  Grayling  data  set, 
environmental  parameters  were  identified  and  incorporated  into  a  linear 
model  that  described  the  variations  in  o®. 

In  conclusion,  we  have  presented  a  small  portion  of  RCS  data  obtained 
during  extensive  measurement  periods  for  two  clutter  sites,  those  at  Gray¬ 
ling,  MI,  and  Yuma,  AZ.  An  empirical  model  has  been  developed  that  suc¬ 
cessfully  identifies  the  primary  environmental  parameters  affecting  the 
clutter  background  for  both  the  Yuma  and  Grayling  environments.  A 
larger  sample  of  cells  would  be  useful  for  future  analyses.  Further  work  is 
needed  to  determine  additional  parameters  responsible  for  some  of  the 
more  subtle  fluctuations  in  the  RCS  measurements,  especially  for  Grayling 
clutter.  This  analysis,  in  combination  with  the  extensive  database  collected 
at  these  field  tests,  will  be  an  important  tool  for  model  validation  and  for 
analysis  of  MMW  system  performance. 
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Acronyms  and  Abbreviations 

A/D  analog-to-digital 

ARL  Army  Research  Laboratory 

ASL  Atmospheric  Sciences  Laboratory 

DAS  data  acquisition  system 

FFT  fast  Fourier  transform 

HH  horizontal  transmit,  horizontal  receive 

HV  horizontal  transmit,  vertical  receive 

I/Q  in-phase/quadrature 

IR  infrared 

LHCP  left-hand  circular  polarization 

LL  left-hand  circular  transmit,  left-hand  circular  receive 

LR  left-hand  circular  transmit,  right-hand  circular  receive 

MMW  millimeter  wave 

PDM  polarization  distortion  matrix 

PRF  pulse  repetition  frequency 

RCS  radar  cross  section 

RHCP  right-hand  circular  polarization 

RL  right-hand  circular  transmit,  left-hand  circular  receive 

RR  right-hand  circular  transmit,  right-hand  circular  receive 

SWOE  Smart  Weapons  Operability  Enhancement 

TIPS  Thermal  Imaging  Processing  System 

V  vertical  polarization 

H  horizontal  polarization 

VH  vertical  transmit,  horizontal  receive 

W  vertical  transmit,  vertical  receive 

WES  Waterways  Experiment  Station 
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Attn46TW/TSWM  B  PameU 
211  W  Eglin  Blvd  Ste  128 
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